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We have found several families of vortex soliton solutions in two-dimensional discrete dissipative 
systems governed by the cubic-quintic complex Ginzburg-Landau equation. There are symmetric 
and asymmetric solutions, and some of them have simultaneously two different topological charges. 
Their regions of existence and stability are determined. Additionally, we have analyzed the relation- 
ship between dissipation and stability for a number of solutions. We have obtained that dissipation 
favours the stability of the solutions. 
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I. INTRODUCTION 

An optical vortex soliton is a self localized nonlinear 
wave, characterized for having a point ('singularity') of 
zero intensity, and with a phase that twists around that 
point, with a total phase accumulation of 2ttS for a closed 
circuit around the singularity [l| . The quantity 5 is an in- 
teger number known as the vorticity or topological charge 
of the solution. Optical vortices can exist in an infinite 
number of ways, as there is no limit to the topological 
charge. This kind of waves looks attractive in future ap- 
plications for encoding and storing information. A spa- 
tial vortex soliton is a specific solution for a (2+1) di- 
mensional nonlinear wave equation [2]. One of the most 
widely used equations of such a type is the nonlinear 
Schrodinger equation (NLSE); it describes wave evolu- 
tion in dispersive/diffractive continuous media with an 
optical Kerr response, i.e. a refractive index that changes 
linearly with the light intensity. When the system under 
consideration has a periodic structure, e.g, a photonic 
crystal fiber, it is necessary to add a periodic transversal 
potential to complete the description, in the NLSE. The 
optical properties of a nonlinear periodic structure can be 
analyzed in the framework of a set of linearly coupled- 
mode equations which, in solid-state physics is called the 
tight-binding approximation, so that the description of 
the system can be understood from a discrete point of 
view. The study of discrete systems has been a hot topic 
in the last years due both to its broad impact in diverse 
branches of science and its potential for technological ap- 
plications Nonlinear optical systems allow us to 
observe several self-localized discrete structures in both 
spatial and temporal domains. 

Unlike conservative systems, self- localized structures 
in systems far from equilibrium, are dynamical solutions 
that exchange energy with an external source (open sys- 
tems). These solutions are called dissipative solitons [7j. 
In Schrodinger models, gain and loss are completely ne- 
glected and the dynamical equilibrium is reached by 
means of a balance between the Kerr effect and dis- 



persion/diffraction. For dissipative systems, there must 
also exist an additional balance between gain and losses, 
turning the equilibrium into a many-sided process 0|. 
The Ginzburg-Landau equation is -somehow- a universal 
model where dissipative solitons are their most interest- 
ing solutions. This model appears in many branches of 
science like, for example, nonlinear optics, Bose-Einstein 
condensates, chemical reactions, super-conductivity and 
many others 0. 

Nonlinear self-localized structures in optical lattices, 
usually referred to as discrete solitons, have been pre- 
dicted and observed for one- and two-dimensional ar- 
rays [TO, EH- The existence of discrete vortex solitons 
in conservative systems have been reported on several 

works EMI. 

For the continuous case, dissipative vor- 
tex soliton families have been found to be stable for 
a wide interval of 5- values [H, EH- Symmetric stable 
vortices have also been predicted in continuous dissipa- 
tive systems with a periodic linear modulation [ItJ . In 
this work, we deal with discrete vortex solitons in dis- 
sipative 2D lattices governed by a discrete version of 
the Ginzburg-Landau equation. We have found differ- 
ent families of these self-localized solutions. We studied 
their stability and found stable vortex families for 5 = 3 
(symmetric) and 5 = 2 (asymmetric) topological charges 
for the same set of equation parameters. In addition, we 
found another symmetric solution in which two topolog- 
ical charges (5 = 2 and 5 = 6) coexist. Finally, we show 
how an increase in dissipation increases the stability re- 
gions for the same "swirl- vortex" soliton analyzed in the 
recent work [l8[ 

The paper is organized as follows. In Section[II]we por- 
tray the model that we are going to use in the rest of the 
paper. Sections IIIII and IIVI describe the new families of 
solutions we obtain, and in Section [V| we compare the re- 
sults of our dissipative model with the conservative cubic 
case (Schrodinger limit). Finally, section PVTI summarizes 
our main results and conclusions. 
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II. MODEL 
A. The cubic quintic Ginzburg-Landau equation 

Beam propagation in 2D dissipative waveguide lattices 
can be modeled by the following equation: 

^m,n + Clpm^n + | ljj m ,n | 2, 0m,n + H VVn,n l^m^n ~ 

+ ^|^m,n|Vm,n • (1) 

Eq.(pQ) represents a physical model for open systems that 
exchange energy with external sources and it is called 
(2+1) D discrete complex cubic-quintic Ginzburg-Landau 
(CQGL) equation. ^m,n is the complex field amplitude at 
the (ra, n) lattice site and ^ m ,n denotes its first derivative 
with respect to the propagation coordinate z. The set 

{m = -M, M} x {n = -TV, TV}, 

defines the array, 2M + 1 and 2N + 1 being the number 
of sites in the horizontal and vertical directions (in all 
our simulations M = N — 8). The tight binding approx- 
imation establishes that the fields propagating in each 
waveguide interact linearly only with nearest-neighbor 
fields through their evanescent tails. This interaction is 
described by the discrete diffraction operator 

where C is a complex parameter. Its real part indicates 
the strength of the coupling between adjacent sites and 
its imaginary part denotes the gain or loss originated by 
this coupling. The nonlinear higher order Kerr term is 
represented by v while e > and /i < are the co- 
efficients for cubic gain and quintic losses, respectively. 
Linear losses are accounted for a negative S. 

In contrast to the conservative discrete nonlinear 
Schrodinger (DNLS) equation, the optical power, defined 
as 

M,N 

Q(*)= E Vv«WV4,„(*) (2) 

m,n=—M,—N 

is not a conserved quantity in the present model. How- 
ever, for a self-localized solution, the power and its evo- 
lution will be the main quantity that we will monitor 
in order to identify different families of stationary and 
stable solutions. 

We look for stationary solutions of Eq. (pQ) of the form 
ipm,n(z) = </>m,n exp(iAz) where (j) m ,n are complex num- 
bers and A is real; also we are interested in that the phase 
of solutions changes azimuthally an integer number (S) 
of 2tt in a closed-circuit. In such a case, the self- localized 
solution is called a discrete vortex soliton [l9| with vortic- 
ity S. By inserting the previous ansatz into model (pQ) we 
obtain the following set of (2M + 1) x (2N+1) algebraic 
coupled complex equations: 

— \<firn,n + ^0m,n + | <t>m,n 1 2 </>ra,n + ^|</>ra,n | 4 </>ra,n = 

iS0m,n + ^|</>m,n| 2 </>m,n + ^|</>m,n | 4 </>m,n • (3) 



We have solved equations (j3j) using a multi-dimensional 
Newton-Raphson iterative algorithm. The method re- 
quires an initial guess, and we have found that usually 
converges rapidly by starting with a high-localized profile 
seed that can be constructed by a procedure similar to 
the one described in [l8| . 

B. Linear stability analysis 

Small perturbations around the stationary solution can 
grow exponentially, leading to the destruction of the vor- 
tex soliton. A linear stability analysis provides us the 
means for establishing which solutions are stable. Let 
us to introduce a small perturbation, 0, to the localized 
stationary solution 

VVn = [<t>m,n + </>m,n 0)] e lXz , (/> m , n G C, (4) 

then, after replacing Eq.(J2]) into Eq.(pQ) and after lineariz- 
ing with respect to we obtain: 

0m, n + Cfpm^n — iS(j) mjn + 

[2(1 - e)|^ m ,„| 2 + 3(i/ - /u)|0 m ,„| 4 - A]^ m ,„ + 

[(1 - e)^ >n + 2(U - »)\4>m,n\ 2 4> 2 m ,n}4>* m ,n = 0- (5) 

The solutions for the above homogeneous linear system 
can be written as 

0m, n (z) ^m,n 

exp [jm,nz] + C^^ n exp [7 m?n ^], (6) 

being C 1 ' 2 integration constants and j m ,n the discrete 
spectrum of the eigensystem associated with ([5]). The 
solutions are unstable if at least one eigenvalue has pos- 
itive real part, that is, if max{Re(7 m?n )} > 0. Here- 
after, we will plot stable (unstable) solutions using solid 
(dashed) lines. 

III. SYMMETRIC AND ASYMMETRIC 
SOLUTIONS 

Eq.([T|) has a five-dimensional parameter space, namely 
C, (5, £, /i, v. In order to look for any stationary solution, 
first, we need to choose a fixed set of values for these 
parameters, and then an initial condition. By starting 
from a guess with eight peaks surrounding the central 
site - the first discrete contour of the lattice around of 
(m, n) = (0, 0)- with a topological charge S = 3 sampled 
on this path, the iterative algorithm rapidly converges 
to a stationary structure with the same features of the 
initial guess. Once we found a stationary solution with 
the desired properties, we use it as initial condition to 
find the corresponding solution for a slightly different set 
of equation parameters. We usually just change one of 
them. Therefore, for the dissipative case we construct 
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Figure 1. (Color online) Q versus e diagram for families A 
and B of discrete dissipative vortex solitons. Inset: Q vs. e 
diagram for family C. (CQGL equation parameters: C — 0.8, 
6 = -0.9, /x = -0.1, v = 0.1). 



families of solutions by fixing four parameters and vary- 
ing the fifth one, usually the cubic gain parameter, e. 
Using this procedure, we have constructed the A fam- 
ily (displayed as the curve Q versus e in FigJT]). We 
started from a highly localized solution and we slowly 
decreased the nonlinear gain, observing that the solution 
became gradually more and more extended as e (and 
Q) decreased. The saddle-node point for this family is 
reached at e w 0.62. 

Figj2] shows the amplitude and phase profiles corre- 
sponding to the solution marked with a green solid circle 
on the A family in FigH]). From the amplitude profile, 
Fig{2^a), we can see how the stationary solution main- 
tains the eight excited peaks of the initial seed. Besides, 
we can see some energy in the tails, i.e. on the second 
discrete contour. On the other hand, the phase profile, 
FigEJb), clearly shows a topological charge S = 3. 
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Figure 2. (Color online) Color map plots for the eight peaks 
stable vortex solution with S = 3, marked with a green circle 
on the A family branch in Fig(T] (a) Amplitude profile, (b) 
Phase profile. 

A similar procedure has been done to construct an- 
other family, labeled B (See FigJT]). This family consists 
of asymmetric stationary solutions characterized for hav- 
ing six peaks displayed on the corners of an elongated 
hexagon in the n-axis direction of the lattice. This spa- 
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Figure 3. (Color online) Color map plots for the six peaks 
stable vortex solution with S = 2 marked with a cyan circle 
on the A family branch in Fig(T] (a) Amplitude profile, (b) 
Phase profile. 



tial configuration possesses a topological charge S = 2. 
Typical amplitude and phase profiles for this kind of solu- 
tion are shown as color maps in Figj3] In the conservative 
case, four peaks structures have been reported to be sta- 
ble [19j for S = 1 and unstable for S = 2; on the other 
hand, for hexagonal lattices, six peak structures are sta- 
ble [20j for S = 2 and unstable for S = 1. In continuous 
systems assymetric four peaks structures has been found 
stable for 5 = 1 HJ. 

The families A and B of stationary vortex solutions 
coexist for the same set of parameters of the discrete 
Ginzburg-Landau equation. Other families of solutions 
exist too for the same set of parameters. The inset shows 
a different family (the C-one), whose Q vs e diagram 
almost coincides with the A-family one, in spite of being 
quite different type of solutions. We will describe these 
C-solutions later in Section IVl 



IV. "TWO CHARGES" VORTEX SOLITON 

Now, we show one example where two topological 
charges coexist in the same solution. Let us start with a 
guess solution consisting of twenty peaks, spatially dis- 
tributed like a rhombus, and with a topological charge 
5 = 2. Using it as the starting point for the Newton- 
Raphson algorithm, we find a stationary solution that 
looks like the one shown in Fig|4^a,b). This solution be- 
longs to the family displayed in Fig |5] labeled with D; it 
was constructed following the same procedure described 
in the previous section. Unlike the previous families, 
the D family does not reach the saddle node point via 
a monotonic decreasing of its power; rather, it passes 
through a minimum value (e ~ 0.64), then, the power 
grows and finally, the saddle node point is reached (See 
inset in Figj5j). 

The solutions of this family present a very interest- 
ing property related to its topological charge. The first 
square contour, Ti, -the innermost discrete square tra- 
jectory on the plane (m, n)- in FigHJb) shows that the 
vorticity has a value 5 = 2. For the second contour T2 
we observe that the topological charge has changed to 
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Figure 4. (Color online) Color map plots for the twenty peaks 
stable two charges (S — 2 and S = 6) vortex solution, marked 
with a green circle on the D family in Fig|5] (a) Amplitude 
profile, (b) Phase profile. 




Figure 5. (Color online) Q versus e diagram for "two charges" 
(S = 2 and S = 6) discrete vortex solitons. Continuous and 
dashed lines correspond to stable and unstable solutions, re- 
spectively. The green circle on the D family corresponds to 
the profiles shown in Fig|4] 



5 = 6. Looking at the remaining contours, we note 
that the topological charge returns to S = 2, so, we 
can talk about a transition of the effective vorticity from 
S = 2—^S = 6—^S = 2^ as we move farther from 
the center. For this reason, we can say that the stable 
solutions of this family possess two topological charges. 
For the sake of clarity, we plot sin(# m>n ) vs ip, the az- 
imuthal angle for the lattice, for the first and second dis- 
crete contours. From Fig|6ja) we can see that the data 
(green points) are perfectly fitted by the sinusoidal func- 
tion (gray line) with two periods (S = 2) along the first 
contour, and for the second contour we have six periods 
(S = 6) as shown in Figj6^b). This is somehow a proof 
of the different topological charges contained in the solu- 
tion, and it also proves that the discrete vortex is a well 
defined structure. 

Fig|5] also shows that the D family has one large sta- 
ble region and another small region, magnified in the 
inset, where the solutions are unstable. These unstable 
structures decay on propagation to another kind of stable 
solutions having less energy, and different amplitude pro- 
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Figure 6. (Color online) sin(# m , n ) versus ip (azimuthal angle 
for the lattice) diagram for the first (a) and second (b) discrete 
contour for the vortex soliton, marked with a green dot on the 
D family in Fig(5] 
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Figure 7. (Color online) Q versus e diagram showing the 
transition from the unstable solution marked with the gray 
circle to the stable solution marked with the green circle; the 
inset shows the power evolution for this transition. 



file with only four peaks. In particular, FigJTl illustrates 
how the unstable solution marked with a gray point (the 
saddle-node point for the D family in the FigH]) decays, 
by means of a radiative process shown in the inset, to 
the stable solution marked with a green circle on the 
E family. The amplitude and phase profiles showed in 
FigsJH^a-b) and FigsIHKc-d) correspond to the unstable 
and stable solutions marked with gray and green circles 
in FigEl 
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Figure 8. (Color online) Color map plots for discrete dissi- 
pative solitons. (a) Amplitude profile and (b) phase profile 
for the unstable twenty peaks vortex solution localized on the 
D family at the gray circle in the FigJT] (c) Amplitude pro- 
file and (d) phase profile for the four peaks soliton solution 
localized on the E family at the green circle in the Fig[71 



V. DISSIPATION AND STABILITY 

In this section, we are interested in analyzing how 
the stability of the solutions is affected when our model 
slowly goes to the Schrodinger limit, i.e when the value of 
the parameters in the CQGL equation (pQ) tends to zero: 
{J, e, /i, v} — >• 0. In particular we will focus on the solu- 
tion marked with a purple circle on the C family in the 
inset of FigJTJ We will compute its stability region when 
the gain, loss and higher order Kerr terms are gradually 
suppressed in the Ginzburg-Landau model. 

These solutions are of the "two-charges" vortex type, 
with charges 5 = 1 and S = — 3 [HI- Moreover, this type 
of solutions (the swirl- vortex soliton) can be understood 
as a bound state of five vortices J2^,j23j. Indeed, we can 
identify a vortex with S = 1 at the origin (O symbol), 
surrounded by four vortex, each with 5 = — 1, whose sin- 
gularities are located at the center of the O symbols on 
the FigEfb). This interpretation agrees with the tran- 
sition of the effective vorticity from S = 1 — 5 = — 3, 
as we move farther from the center. The amplitude and 
phase profiles for this solution are displayed in FigJH^a- 
b). It should be noted that it has the same power value 
that the solution marked with the green circle on the A 
family. (In fact, we plot its corresponding family in the 
inset because both families have almost identical Q ver- 
sus e diagrams). We can see that both amplitude profiles 
are very different; even though each one has eight princi- 
pal excited sites, their spatial distributions are dissimilar. 
In addition, their phase profiles are completely different. 
While the solutions on the A family have a well defined 
unique topological charge 5 = 3, the solutions on the C 
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Figure 9. (Color online) Color map plots for the eight peaks 
stable two charges (S = 1 and S = — 3) vortex solution lo- 
calized on the C family at the purple circle in the Fig[T] (a) 
Amplitude profile, (b) Phase profile. 




Figure 10. (Color online) sin(# m ,n) versus ip (azimuthal angle 
for the lattice) diagram for the first (a) and second (b) discrete 
contour for the swirl- vortex soliton, marked with a purple dot 
on the C family in Fig[T] 



family have two charges (5 = 1 and 5 = — 3) simulta- 
neously, as mentioned before. Again, we plot sin(# m)n ) 
vs cp for the first (Ti) and the second (I^) discrete con- 
tour. From Fig fTUT a) we can see one period (5 = 1) 
for the sinusoidal function (gray line) along the first con- 
tour, and for the second contour we have three periods 
(5 = —3) as shown in FigfTUTb). Unlike the conservative 
cubic case (NLSE), in the dissipative model the prop- 
agation constant A is not an arbitrary parameter that 
can be chosen at will. It is fixed by the rest of the CQGL 
equation parameters. By changing them, the value of the 
propagation constant also changes. As in other nonlinear 
problems [2J], we can think of the dissipative terms as 
determinant to select one of the infinite solutions of the 
associated conservative problem. With this in mind we 
will find out the stability regions in terms of the propa- 
gation constant so we can compare with the Schrodinger 
limit. 
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For the sake of comparison we construct the Q versus 
A diagram shown in FigJTTJ Here, we have fixed J, \i and 
v parameters and we only move through the e parameter 
(nonlinear gain). In this way, we obtain a solution and 
its corresponding propagation constant for each value of 
e. Then, we proceed varying the rest of the parameters 
slightly, and construct a new curve, taking the solutions 
of the previous curve as initial conditions in our mul- 
tidimensional Newton- Raphson scheme. In the inset of 
FigJTT]we show the corresponding A vs e diagram. 
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Figure 11. (Color online) Q versus A(e) diagram for several 
sets of parameters specified in Table I, of two charges swirl- 
vortex solitons. Inset shows A vs £. 

With the previous scheme we can find a large number 
of curves, but for the sake of clarity, we only show three 
of them; they are located between the conservative cubic 
case (black branch) and the C curve (gray branch). We 
can read from Table U the CQGL equation parameters 
corresponding to the curves displayed in FigJTTJ These 
five branches belong to the same family conformed by 
vortex solutions with amplitude and phase profiles such 
as those showed in Figj9] 

We have done the standard linear stability analysis, 
described in section II, for each one of them. The part of 
the curves that correspond to stable solutions are shown 
with continuous lines while the dashed lines correspond 
to unstable solutions. 



Taking the above into account, we can clearly establish 
that if the dissipation is attenuated the stability regions 
for the soliton solutions are reduced. [24| Indeed, we can 
see here a wide difference between the stability regions 
for the Schrodinger limit and the C branch. The first 
one only has stable solutions for propagation constant 
values far away from the linear band, the last one has 
stable solutions for propagation constant values closer to 
the linear band. 

Table I. CQGL equation parameters 
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VI. SUMMARY AND CONCLUSIONS 

In conclusion, we have found discrete vortex solitons 
(symmetric and asymmetric) with higher-order vorticity 
in dissipative 2D-lattices and studied its stability. In par- 
ticular, we have also shown in detail a solution that con- 
tains two topological charges. Finally, we analyzed the 
stability of the solutions when dissipation in the system 
is decreased, observing that the stability regions shrink. 
A comparison with the conservative cubic case is done, 
showing that dissipation serves to provide stability to 
otherwise unstable conservative solutions. 
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